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In view of recent developments in the investigation on cuprate high-T c superconductors and the spin- 
Peierls compound CuGeOa, we study the effect of dilute impurity doping on the spin-Peierls state in 
quasi-one dimensional systems. We identify a common origin for the emergence of antiferromagnetic 
order upon the introduction of static vacancies, and superconductivity for mobile holes. 



The evolution of a disordered spin-gapped state (here- 
after referred to as a spin liquid (SL)) into a superconduc- 
tor as observed in the underdoped regime of the cuprate 
oxide compounds, continues to pose a major mystery. 
Here we address this problem from a new perspective; 
the comparison of the effect of static and mobile vacan- 
cies. 

Doping the CuC>2 planes with static nonmagnetic im- 
purities provides valuable information on the SLs. While 
monovalent Li can lead to additional complications [Q, 
Cu— »Zn substitution which introduces vacancies without 
excess holes serves as an ideal probe: enhancement of 
antiferromagnetic fluctuatuation accompanying the nu- 
cleation of local s—1/2 moments is observed ap- 
parently indicating that the state is a spin singlet with 
confined spinon excitations. Disorder-induced antiferro- 
magnetism (AF) is a property also shared among several 
quasi-ld spin-singlet systems, most notably in the spin- 
Peierls compound CuGeC>3 jl , intensively studied in the 
past several years ||. Such analogy has attracted much 
attention; the effect of disorder on the density of states 
of the staggered flux state has been discussed B in this 
light, as well as a possible explanation |7| of the extreme 
sensitivity of the 40 meV magnetic resonance peak to Zn 
impurities. However the precise relation, if any, between 
such responses of SL against static impurities, and the 
superconducting instability observed in the cuprates (or 
spin-ladders) in the presence of mobile holes remains un- 
clear. This is the issue we wish to take up below. 

In this article we have chosen to re-examine the sim- 
plest confining SL, the quasi-ld spin-Peierls (SP) sys- 
tem within the above scope. The primary aim is to in- 
voke nonperturbative methods which takes advantage of 
the one-dimensionality, and to extract a physical picture 
which may well be generic to a wider family of confin- 
ing SLs. Additional motivations come from experiments, 
however; while carrier doping so far has not been real- 
ized in CuGeOa, a possible proximity/coexistence of d- 
wave superconductivity with a SP-like (or bond-centered 
density wave) state in underdoped cuprates Q has been 
inferred by very recent neutron scattering data of longi- 
tudinal optical phonon dispersions in Lai.ssSro.isCuO^ 
and YBa 2 Cu 3 6 .95 ||. 

Let us summarize in physical terms what will follow. 



Theories on impurity effects in spin chains have often 
focused on either the formation of local spin moments 
(using AB), or on quantum interferences among a pre- 
assigned array of antifcrromagnetically aligned spins in 
the presence of a few vacant sites (employing semiclassi- 
cal methods). The former cannot account for directional 
fluctuations of the spins, while the latter may overem- 
phasize the role of local spin moments by assuming them 
even in a spin singlet state. Our remedy is to devise a 
version of nonabelian bosonization (NAB) which simul- 
taneously resolves both aspects. We derive an effective 
action reproducing the known AB theory for SP systems 
|iof , with two additional terms, each related to the di- 
rectional fluctuations of the AF order parameter and the 
Berry phases. The contribution of these new terms are 
appreciable only near an impurity- induced spin moment. 
This local weight transfer together with the Berry phase 
effect (which keeps the spin moments in registry with the 
underlying AF pattern) is what gives rise to AF order in 
the static impurity case. An important aspect here is the 
length (or energy) scales involved. Adiabatic adjustment 
of the spins to the charge deficiency requires the healing 
length of the charge to be sufficiently shorter than that 
for the spins. Meanwhile when mobile holes are doped, 
the same basic conditions are seen to enhance supercon- 
ducting instability, enabling us to make contact with a 
pairing picture proposed in several earlier works. In this 
case, the local AF enviroment provided by the spectral 
weight transfer enhances intrasublattice hopping as well 
as mediates intersublattice attraction. The adiabadicity 
condition ensures slow fluctuation of the AF enviroment, 
necesarry for coherent motion of the holes. 

To model a SP system, we incorporate the one-band 
Peierls-Hubbard (PH) model at half-filling, 

H = J^it ~ (-1)^(4^ + h.c.) + C/^n lT n lX , (1) 

icr i 

with U > 0. Interchain coupling (which fixes the pref- 
erence of the dimer pattern and hence confines disorder- 
induced spinons) will be assumed. We start by a semi- 
classical description of the bulk state; treating the U- 
term as a commensurate spin density wave (cSDW) sad- 
dle point solution (p =< c) ^c >, with 0i = (— l) l mfii 
(|n,:| = 1, m ~ jje^TT 1 [pi), we get a 4-component dirac- 
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fermion type Hamiltonian density E3] 
H F = [R\L^] 
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(9±) and spin (<fr±) phase fields. We attempt a mod- 
ification by the simple replacement a — > Q = ft ■ a. 
Constructing fermi-bilinears according to this rule, one 
readily sees that it corresponds to parametrising the k=l 
Wess-Zumino-Witten (WZW) field g GSU(2) appearing 



where the right (R) and left (L) movers each carry a spin in the NAB rules [|13j as g = e 
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(The dual fields 



index, Q = n-a, A = y / («^»)2 + ( 2 5t) 2 , and tan(^) = 

i^j^- To see the physics embodied in the mass term 
(off-diagonal elements) assume temporarily that ft = z. 
The effective spin-dependent potential energy V a (x) = 
e t2kFX V a (2k F ) + e- l2kFX V a {-2k F ) experienced by each 
spin component a = ±1 can be read off using W ff-diag — 
Z a RiLM2k F )+LlR a VA-2k F ): 

vM _f 2A sm(2k F x + (§ + &)) (a = +l). . . 
a{) \ 2A sin(2A; Fa: -(f + f )) (a = -1). W 

The minima of the potential are located at x m i„ = Xa 
(a: lattice constant), where X = (2n — 1) — j^n € Z 
for a = +1 and X = 2n + |j,n € Z for cr = -1, which 
invites the following interpretation. When 4> = 0, odd 
sites (even sites) are occupied by down-spins (up-spins) 
(the cSDW theory). Turning on electron-lattice coupling 
(4>a =/= 0) shifts the position of these down-spins (up- 
spins) to the right (the left) resulting in a regular array 
of strong (odd-even) and weak (even-odd) bonds. The 
antiparallel spin pairs on the strong bonds form bond- 
centered density waves. Returning to the general case, 
this picture remains locally valid with the replacement 
z —t n provided ft fluctuates on a scale longer than a. 
A low energy theory is obtained by treating the corre- 
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'7o, 7 5 = «7o7ii and e lQ * 7 
1 + i sin 4r<3 ® 7 5 ) in a derivative expansion of ft. Per- 
turbative terms enter only beyond quadratic order [fl3|| , 
and the resulting action - coming from the anomaly of 
the SU(2) current = ^j m j 5 ^ - is the 0(3) nonlinear 
sigma (NLer) model with vacuum angle 9 — n— (f>o — sin^o 
[|l2],0, which signals a spin-gap for O ^ 0. 

Having discussed the bulk system at 1/2-filling, we now 
introduce a dilute density of nonmagnetic impurities (va- 
cancies). To this end, keeping </>o constant (justified in 
the presence of interchain coupling), we bosonize. In do- 
ing so we must go beyond usual AB with a fixed spin 
quantization axis in order to retain the spin directional 
degree of freedom ft. We describe in some detail how this 
can be done. The AB scheme derives from the fermion- 
boson correspondence 
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(in conventions of ref. |15|), and translates fermionic op- 
erators into the language of a set of conjugate charge 



and 0- are gauge degrees of freedom which do not arise 
unless considering chiral currents or Cooper channels). 
Note that we differ from the usual way JO] of relating 
abelian phase fields to the WZW model. To give firmer 
grounds to this identification observe that eq.(|J) and its 
nonabelian generalization may be viewed as a family of 
chiral transformations acting on a bosonic vacuum free of 
charge or spin solitons. The spin part of the free fermion 
theory therefore has an induced SU(2) connection and is 
consistently evaluated as 



= / VnD4>+e 
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where U§ = e~^ +< ^ 75 . The bosonization dictionary re- 
mains unaltered for the charge sector, while the spin part 
receives corrections related to the fluctuations of ft. For 
instance the k = 2k F component of the spin operator 
reads Sk=2k F oc sin(2k F x + 9 + ) sin<^>_|_n, whereas the uni- 
form component is 
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in which ^additional = ^ cos (f> + sin (j) + d x n- ^ sin 2 <p + n x 
d x ft. Canonical quantization of ft and 0+ yields the 
correct (Kac-Moody) algebra for the currents Jr and 
Jl- The bosonized Lagrangian for C F is £ = C{9 + ) + 
£wzw(g)| s = CX p(-i0 + Q) + C m ass(0+, <f>+,ft), where the free 
fermion parts for the charge and spin are respectively 
C e+ = ^9+f and 
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£wzw(5)l S =c X p(- l( /. + Q) = —(d M <f> + ) +^ sin <M^«) 
+^(20 + - sin(20 + ))g TX , (7) 

with q TX = j^ft ■ d T ft x d x ft, the topological charge den- 
sity of instantons. The last two terms in eq.(0) are new; 
they are terms that disappear on taking n=const., repro- 
ducing the AB expression for a free Tomonaga-Luttinger 
liquid. On the other hand, a bulk spin gap will fix </> + , 
and (for cf> + ^ 0) these terms will yield an 0(3) NLct 



model with a (9-term. The mean value of 



is deter- 



mined from the interaction C mass , but the latter needs 
be handled with some care. Straightforward bosonization 
gives 
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where a is a short distance cutoff. Being a relevant 
term, this locks 0+ in the bulk problem to the value 
+ = § — ^f- (We assume that sin6* + attains a mean 
value, i.e. the Umklapp term present at half-filling opens 
up the charge gap.) Plugging this into £ wzw , we re- 
cover the NLcr model with the 0-angle previously men- 
tioned. But to see the full correspondence with the 
AB result, we should go beyond this semiclassical ap- 
proximation; noting that the portion of the mass term 
Ao cos ^{RaQapLp + h.c.) had originated (prior to the 
decoupling) from a backscattering process, we should cor- 
rect eq.(@) into the form 
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where we now have complete agreement with the well- 
known AB result for the spin-Peierls system |to| , supple- 
mented with the second and third terms of eq. @) . The 
D-term is marginally irrelevant, and the effective value of 
9 is now governed by the first term of eq.(||). For 0o ^ 
this gives 0+ — and hence 9 — 0, in which case the 
magnitude of the staggered spin mft is quenched, mak- 
ing the second term in eq.(Q) ineffective. (The formula 
for Sk=2k F infers that 9 = and tt each corresponds to 
a spin-singlet and a Neel state.) The case 0o = (no 
dimerization) is special; only the D-term is present and 
the effective 0-angle is undetermined, indicating a dy- 
namically induced axial U(l) symmetry. This suggests 
a physical picture of fluctuating dimers reminiscent of 
a long-ranged resonating- valence-bond state j[o|. The 
above arguments expose an intimite and rather unex- 
pected relation between the 0-angle and the spin phase 
fields of AB, which only becomes evident in the present 
"rotating frame" . In the remaining part we seek its con- 
sequences, concentrating on the case 0o 7^ 0. 

We are now ready to discuss vacancies. From the for- 
mula for the charge density p = ~d x 9+, this is repre- 
sented by a 7r-kink of 9 + . From either eq.(||) or eq.(||), 
this is seen to invert the sign of the potential energy, 
which must be compensated by a 7r-kink of 0+ . The lat- 
ter corresponds, according to eq.(|^) to the liberation of 
a spin 1/2 degree of freedom in the background of the 
singlet state. {Note that this argument does not apply 
in the absence of the spin gap, i.e. for 0o = 0.) ft is 
possible to show, that along the lines of ref jl6| the spec- 
tral weight Imx(k,ui) for fixed 0+ can be estimated as 

2 i 

~ , cos + for the gapped part, while another contri- 

bution ~ | sin0 + /(fc— 7r)| represents a spin wave-like part. 
Hence the 7r-kink of + should indeed cause a transfer of 
spectral weight into subgap states, which is a characteris- 
tic feature of the present system. Physically, a vacancy at 
site x = Xi should release a spin of S ~ |(— l) Xi / a ft(Xi), 
which fixes the sign of the kink of the + field to be 

W+ = TEx«(-l)* /a 9(s--Xi) = {-l) X ' /a 0+. (6 is 
the step function.) This should be true under the adi- 



abadicity condition £ c <C £ s where £ c (£ s ) is the charge 
(spin) correlation length. The continuum limit however 
does not distinguish which of the two sublattices a given 
point x = Xi belongs to. Such lattice effects can be par- 
ticularly important when dealing with Berry phases [JL7| . 
To continue working in the continuum, we are thus lead 
to introduce two charge phase fields 9+ and 6+, one for 
each sublattice. This leads to a simple expression for the 
deviation 5(j>+ of 0+ from the bulk value + , 
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which is the principal equation of this article. Now let us 
see how this affects the topological term, C top = i [20+ — 
sin(2(j) + )]q TX . Again using £ c -C £ s , the effect of <50+ on 
the term — i sin(20 + )c/ T:E cancels out on average, and 



Ctop = i[2<p+ - sin 2cj> + ]q TX + 2i[9% - 9f}q T 
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provided the average seperation of vacancies I > £ s . 
Next, we note jl|] that formally q TX = -^d x Ao, where 
Aq(t, x) = d T ft ■ a(n(r, x)), and the monopole vector po- 
tential a satisfies x a = ft. Integrating by parts, the 



second term in eq. (Jill) becomes C' t 
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d T ft. For static vacancies, this yields the action S' top = 

|Ex 4 (- 1 ) A:</ M^(' r >*i)]) where W [«( T )] = JdrAa is 
the solid angle subtended by n(r) in the course of its 
evolution. These terms are the Berry phases of spin 1/2 
objects induced by vacancies. Together with the bulk 
contributions consisting of the NLcr model and the first 
term of eq.(ll) (we use for the potential eq.(^)), this is 
essentially the action derived in ref. |nj for the doping 
of a spin ladder. Following similar arguments, we arrive 
at the final action for the induced spins {npQ)}, 

S eff [{n j }]=J2k-l) Xj/ Mnj(T)} 



J dTj eff e^- x ^/^n 3 (T) ■ n 3+1 (r), (12) 



where n, = ft(Xj ) and J e f / = l/£ • sin 2 ^ . Absorbing 
the signs into the spins Nj = (— l) x ^ a rij, this becomes a 
random exchange Heisenberg model, with diverging spin 
correlation and staggered susceptibility at T=0 |l^] . We 
expect that the essential physics of the disorder-induced 
AF observed in CuGeC>3 is captured within this model. 

Turning to the case of mobile vacancies (holes), the 
part of the action involving 9+ and re 
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This coincides with the action proposed by Shankar |2Cj ] 
on semiphenomenological grounds for hole motions in an 
antiferromagnetic background. We have arrived at this 
form from an electron system containing both spin and 
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charge sectors. Hereon we may basically adapt the argu- 
ments of ref. pel . Refermionizing £(0+,0+) we see that 
it is equivalent to two massless fermions 



the conference RPMBT11. 



£hoie = + iA)ip A + iPb{0 - iA)i>i 



(14) 



coupled to the gauge fields = n ■ d^a each describing 
intrasublattice (next nearest neighbor) hopping of the 
holes [^,^2|. Because the fermions ipA and ipB have op- 
posite gauge charges, their is an attractive interaction. 
The spin singlet superconducting susceptibility is the cor- 
relation fuction of (in terms of the original electrons) 
l&iW L-*(*) ~ ei ( - e +- e ^e^ e - +e ^e^+. The non- 
trivial combination here is 9+ — 9+, but Gauss law con- 
straints can be incorporated |20|] to show it is massive. 
Then the susceptibility obeys a power law with exponent 
-1, and should become an Emery-Luther superconduc- 
tor when including interactions among neighboring sites, 
which seems to be consistent with available numerical re- 
sults on dimerized t-J models p3|| . This would suggest 
the possibility of the coexistence of the spin-gapped state 
(i.e. spin-Peierls state) with superconductivity. Super- 
conductivity in quasi-ld spin-gapped systems may be- 
come relevant in view of the recent advances in hole- 
injection via field-effect transistors plj . Finally we note 
certain differences from ref [^0| ; first the picture relies on 
a confining SL, therefore breaking down at the Heisen- 
berg point cj>o — 0, i.e. a spin gap is required. A second 
feature is the difference in the vacua structure; we had 
effectively 9 = irrespective of the value of (f>o ^ 0, and 
hence zero weight for the NLtr model part, so the collapse 
of a 6*-vacua structure with hole doping [^o| is not seen, 
which marks a departure from semiclassical methods. 

The crux of our argument is that the enhancement of 
superconducting susceptibility, related to the pairing sce- 
narios based on t'-J type interactions (22|, has emerged 
here from the same origin as the disorder-induced AF, 
which is the coupling of the charge density fluctuation 
to the spin gauge fluctuation, represented compactly in 
eq. ([h]) . This suggests that the validity of such pictures 
in a particular spin-gapped system can be inferred from 
testing its response to static nonmagnetic impurities. In 
this respect we mention recent hole-doping experiments 
in the Haldane-gap material Y2BaNiC>5 p4}| , a system 
where static impurities do not induce AF. The authors 
find no enhancement of conductivity and ascribe it to the 
comparable magnitude of the charge and spin gaps. 

Since our framework can be readily applied to charge 
stripes modeled as arrays of Id electron gases coupled to 
spin-gapped chains ]2q ], similar treatments should pro- 
vide useful insight. Such work is now in progress. 
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